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ABSTRACT 



^ . We extend our earlier work on static color singlet correlators in high T 

QCD (DeTar correlators) to baryonic and gluonic sources, and estimate the 
corresponding screening masses using the dimensionally reduced theory. We 
discuss spin and polarization dependence of meson and baryon masses in the 

OO ■ T — i> oo limit, and possible nonperturbative effects at non-asymptotic temper- 

atures. 
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1. Introduction 

Recent lattice measurements of correlation functions of the type 

C^{x3) = {f^dT f£r.Ur,Xs,T) T t/r'J„(0, 0, r')) (1.1) 

JO J Jo 

(where /5 = 1/T, r = (xi,X2) and Jq, a color-singlet operator)[| show that the 
hadronic screening lengths above the QCD finite temperature (chiral) phase 
transition point fall into chiral multiplets. The screening masses are approxi- 
mately 27rT for mesons and SvrT for baryons, with the notable exceptions being 
the pion and its scalar partner, which are considerably lower |l|. It has also 
been known for some time that bulk thermodynamic quantities, like energy 
and entropy density, are roughly as expected from a gas of non-interacting 
quarks and gluons [§. The exception here is the pressure, which is consid- 
erably lower than expected from the free gas picture. Finally, the fermionic 
susceptibility, d^ In Z/dm? varies rapidly across the transition region from be- 
ing small in the hadron phase, to being large in the quark phase 0. All these 
results are suggestive of a high temperature phase of weakly interacting (Debye 
screened) quarks and gluons. 

On the other hand, the Coulomb gauge lattice calculations of the spatial 
structure of the mesonic correlators via 

Ca(^,a;3) = ( / dr j dVi dVs ^'"^(fi, 0, r)r«^(f2, 0, r) 

r dr' ij\R, X3, T')T^ij{R - f, X3, t')) (1.2) 

JO 

show strong evidence for correlations in the transverse direction [^ . The pres- 
ence of such correlations is not unexpected given that the the spatial Wilson 
loops obey an area law at all temperatures 1§],|§[,[0. 

What causes the correlations in the spatial directions, and to what extent 
they are important for our understanding of the high temperature phase, is not 
clear. What is clear, however, is that any description of the high temperature 
phase has to account for the above results. In a recent paper [^], hereafter 
referred to as I, two of us suggested that at very high temperature, the screen- 
ing lengths and "wave functions" discussed above, can be understood from an 
analysis of the static part of the gluon field together with the lowest energy 
quark modes. This " dimensionally reduced" QCD is a YM-Higgs model with 



^ For baryonic sources the average over imaginary time, r, is weightened by cos loqT where 
L0{) is the lowest Matsubara frequency. 



heavy {i.e. M = ttT) quarks, which is beheved to be confining.^ The screening 
masses and wave functions correspond to masses and wave functions of heavy 
quark states in the dimensionally reduced theory, and can be calculated as 
Coulomb bound states using usual charmonium-type methods. The masses, 
rria, in the mesonic correlators are thus naturally rria ~ 2M = 2itT, and the 
wave functions show exponential fall-off at large distances, both in qualitative 
agreement with the numerical calculations. It has also been shown, that in 
the case of 2+1 dimensional QCD at high T, the Coulomb bound state pic- 
ture leads to the correct description of the quark susceptibilities at asymptotic 
temperatures P|,p!0[]. 

In this paper we extend our analysis in I to (color singlet) baryonic and 
gluonic sources and give a more detailed discussion of the spin and polarization 
dependence of the screening masses. For the baryons the calculations are 
very similar to the mesonic case and are on the same level of rigour. In 
the case of gluonic operators, our methods are more questionable since the 
size of the resulting bound state is larger than the Debye screening length, 
and we cannot rule out the possibility that the infrared divergences in the 
magnetic sector will change our predictions considerably. In section 2 we 
recapitulate the arguments that lead to the dimensional reduction scheme 
using a formalism that is somewhat different from the one used in I. We also 
discuss the polarization dependence of the various mesonic correlators. Our 
results for baryonic and gluonic correlators can be found in section 3 and 4 
respectively. In the latter we also consider correlations between Wilson lines. 
In section 5, we discuss non-asymptotic effects that might explain the difference 
between the lattice data and our asymptotic prediction for the fine-structure 
of the screening mass spectrum. Section 6, finally, contains some summary 
comments and discussion. Some details about the calculations can be found 
in Appendix, A, B an C. 

2. Dimensional Reduction and DeTar Correlators 

We now briefly summarize the dimensional reduction scheme for the fermionic 
part of the high T QCD lagrangian. The starting point is the Euclidian action 
for free fermions, 

C^ = itP^-%D^ + m)tij , (2.1) 

where D^ = d^ — igA^., and where the (Hermitian) Euclidian gamma matrices 



^ This is based on the already mentioned fact that space-hke Wilson loops obey an area 
law at all temperatures. Recent work by Karkkainen eta/.lq] shows that the string tension 
cr(T) ^ T^ at high temperature. 



satisfy 

{7^,7.} = 25^, (2.2) 

Now, retain only the modes with energy uj = ±7r//5 = ±7rT and define (f)± by 



where U± = e^' 2"^ with cot i!) = m/nT, and where V = e*""''^'^/^ with a^j^ = 
|[7/i,7i/]5 rotates an angle 7r/2 in the (/ii/) plane. In particular \^"'^73\^ = 74 
and V^^4V = —73. Expressed in these variables the dimensionally reduced 
Lagrangian becomes 

£3 = Y.-z^l[%03-tgAs)+^,{di-tgA,)-M (2.4) 

+ ig A4{cos 1!}% — ai sin ^)](l>a 

where M^ = {ttT)^ + 771"^. The last step is now to rotate to a fictitious 2+1 
dimensional Minkowski space (from now on greek indices run from to 2 and 
roman from 1 to 2) by = cj), (p^ = icj), {xi,Xs) = {x\ix^), (71,73,74) = 
(-i7j, -^73, -70) and (1^,^3,14) = VT{Ai,-iAo,H), where H is the Higgs 
field, gs = VTg the 3d gauge coupling constant and where we use the metric 
(H ) in Minkowski space. The resulting Lagrangian readsQ 

/:2+i =T.^a ^7''('9m - tgsA^) -M + gsiasin^ - cos ^7')^] <Pa • (2.5) 

a=± 

In the dimensionally reduced theory the quarks are heavy with mass M, the ad- 



joint Higgs scalar has the electric mass niE = j4:/3gT (for two light fermions) 
and the gluons are expected to acquire a magnetic mass ~ g^T. The quark and 
Higgs masses are easily understood in perturbation theory as a consequence 
of the absence of fermionic zero modes and Debye screening respectively. In 
contrast, the magnetic mass is not calculable in perturbation theory and the 
predictions made on the basis of the dimensionally reduced theory are only 
reliable if they are insensitive to distance scales > 1/g'^T. A more detailed 
discussion of these points can be found in I. 

■^In I we used a different reduction sclieme based on two-spinors. Tlie lagrangian given 
there was not correct due to a subtlety in one of the transformations. However, the conclu- 
sions of the paper are not affected. We give the correct two-spinor lagrangian in Appendix 
A. 



Note that our conventions are slightly unusual in that we use 4 dimensional 
matrices to represent the 2+1 dimensional Clifford algebra. In Appendix A 
we give a formulation in terms of two different two-spinors that however have 
non-diagonal mass terms. For practical calculation the Lagrangian (2.5) is 
more useful. 

For m = the original Lagrangian (2.1) is invariant under the chiral trans- 
formation ip -^ e'^'^^'^ip, ip'^ -^ %l)'^e^"-^^ where 75 = 75 = —71727374- After the 
variable change, this corresponds to the Lagrangian (2.5) being invariant un- 
der 0-1- -^ e^*°'''^''''5 0±. Also note that in spite of the mass term £2-1-1 must be 
invariant under parity transformations since the original £4 is. It is easy to 
verify that the relevant transformation is Xi — > — Xi, (j)± — > —^75710^:. 

Under the variable change (2.3) the expressions for the sources are also 
changed and we give a translation table for the most important currents in 
the chiral limit tti = 0. Others are easily arrived at using (2.3): 



i,^ - 


^ +#±730± 


^Hs^ - 


-^ ^0±750± 


ii^4) - 


-^ +0±0± 


^Hs^ - 


-^ -#±7o0± 


^Hi^ - 


^ 0±7j0± 



(2.6) 



As explained in I, the screening lengths and wave functions for mesonic 
and baryonic operators, can be identified with the masses and wave functions 
of the non-relativistic bound states corresponding to (2.5). The relevant non- 
relativistic Hamiltonian is obtained by Breit reduction. Keeping only the 
leading Coulomb interaction we have 



* = i:|7+i:|^(i + >n(Mh-r,|)) (2.7) 

where e^ = g\Cp (quark-antiquark) and e^ = g\Cpl1 (quark-quark) with Cp 
the Casimir operator in the fundamental representation. To arrive at this 
expression one has to cancel the infrared singularities between self-energy and 
exchange contributions as explained in I. 

We now give a formula for the spin-spin interaction in non-relativistic 
bound states by performing a Breit reduction of the Hamiltonian correspond- 
ing to (2.5). The easiest way to proceed is to notice that the expressions for the 
currents are exactly as in the standard (3+1) dimensional case, except that the 
momenta in the 3-direction are identically zero. Since, in the nonrelativistic 



approximation, the spin-dependent part of the current is x^(p')^o^x {f^~p)x{p)j 
the Pauh interaction only depend on the spin component in the 3-direction. 
A simple calculation yields, 

^s. = E 717^(1 + ^SlSi) 5\n - r,) (2.8) 



i<j 



4M2 



Hss is subleading in the temperature and, as usual, it will only be considered 
as a perturbation. In the case of mesons (2.8) simplifies to 

Hss = ^,Sl6\n-r,) , (2.9) 

where 5*3 is the total spin of the meson in the 3-direction. As mentioned ear- 
lier, this expression differs from the one derived in I. In the final result there 
is, however, only one difference, namely that the different polarization compo- 
nents of the rho acquire different screening masses. Of the four components of 
the interpolating massive vector field p^ = ipj^iip, the components 1, 2 and 4 



are measured on the lattice [|Tl].0 Using the translation table (2.6) we see that 



these correspond to (frfi(j) and (j)(j). The first has spin-projection ±1, while the 
second is not an eigenstate of 5*1, but a mixture of 1 and 0. The component 
Ps, which has not been measured on the lattice, corresponds to the current 
07o0, and is a 5*1 = eigenstate. 

Using (2.9) we predict 171.^ — ^ttT < rup^ < rup. and | 

The lattice measurements at temperatures not much above Tc do not exhibit 
this pattern but rather |m^ ~ nip- ~ 2'kT < nip^. One should remember, 
however, that our predictions are good at asymptotically high temperatures, 
while the lattice calculations are performed not far above the transition tem- 
perature. The consequences of this will be discussed further in section 6. 

At this point it is also pertinent to make a few comments on the flavour 
assignments in the lattice simulations of the mesonic screening lengths that we 
have reffered to. These were carried out using staggered fermions, and in this 
scheme the space-time and flavour symmetries are intermingled on the lattice 



■* In ||ll| Pi and p^ are referred to as vtl and utO respectively, where the and 1 denote 
the " hehcity" which is the component of the angular momentum about the 3-direction. We 
do not agree with this helicity assignment for /94, as discussed in the text. 

^ Similar results have also been obtained by Brown eta/.using somewhat different 
methods, ffll 



and are retained only in the continuum limit. The pion and sigma quoted 
originally by Detar and Kogut are related to an exact ^(1) ® U{1) symmetry 
of the free lattice fermion theory. The generators corresponding to sigma and 
pi are 1 ® 1 and 75 ® T5 respectively, and the near degeneracy of their screening 
masses above T^ is interpreted as a restoration of the symmetry. However this 
may be subject to some doubt, since, as first pointed out by Shuryak |jl3[, only 
the connected part in Fig. 2a was simulated on the lattice. The disconnected 
part in Fig. 2b remains to be calculated. 

Since there is no lattice measurements of any other member in the pseu- 
doscalar multiplet we still do not know whether the whole continuum 5'[/(4)ij{8) 
SU{4:)l is restored, or only the above mentioned subgroup. There is also no 
information on whether or not the axial Ua{^) symmetry (not to be confused 
with the lattice symmetry mentioned above) is restored. This question, which 
is directly related to the strength of the anomaly at finite temperature is impor- 
tant in order to correctly identify the symmetry group of the high temperature 
phase, and is thus crucial for arguments of the type recently advanced by Ra- 
jagopal and Wilczek [|1^]. Clearly it would be very useful to perform a direct 
finite temperature lattice calculation of the f7A(l)-types of order parameters 
such as 



det(v^+^(l±75)V^) 



as originally suggested by t'Hooft |15 



3. Baryonic Correlators 

The DeTar correlator in the baryonic channels can be calculated completely 
analogously to the mesonic ones. The currents appropriate for the nucleon and 
the delta isobar, omitting color, flavour and space indices, are 

J^ = {^'ilj^C%^)ij + K{tp''Cij)%i^ (3.1) 

J^ = {^''C%^P)^ (3.2) 

where k is an arbitrary mixing factor and C = 7274 the charge conjugation 
matrix. The existence of the two nucleon currents (hence k) was first noted 
by loffe [|l^. To use the dimensional reduction method we should transcribe 
(3.1) and (3.2) in terms of the fields 0±, as shown in Appendix B. Using these 
currents, the asymptotic form of the correlator (1.1) gives the screening mass 
rUa in the baryonic channels. To leading order rria = 3M + E^ is just the 
mass of three heavy quarks in a color singlet configuration, interacting via 
Coulomb forces. In the ground state, the three quarks are in a symmetric 



spatial configuration, in which case (2.8) simphfies to 

Hss = ^ {S! + I) 6\n,) , (3.3) 

where S3 again is the total spin along the 3-direction and ri2 the relative two 
body separation in the three body system. This is to be compared with the 
spin-spin interaction (2.9) in the mesonic channel. 

The ground state baryon wavefunctions are antisymmetric in color and 
symmetric in spin, flavour and space, and are solution to the Faddeev equa- 
tions. For simplicity, however, we will only make variational estimates. In 
the center of mass frame we use the following Gaussian variational ansatze for 
the ground state wavefunction (we also give results for high T QCD in 2-1-1 
dimensions^) , 

gCD2+i: ^(p,r^) = W— e-'^(^+'') , (3.4) 

V TT 

QCD3+1: ^(p, r/) = — e-"('''+''') , (3.5) 

TT 

where p = 4s (^1 — "^2) and f/ = w | [fx -|- r2 — Ir^^ . Minimizing the expectation 
value of (2.7) in (3.4) and (3.5) with respect to a yields for QCD2+1, a = 
I f 9eM \ 3 ^^^ £^^ QCD3+1, a = ^j^- The corresponding baryonic screening 
lengths are, 

9 / 36^ \ ^ 
QCD2+1: ^" = 3M+-(^j , (3.6) 

QCD,^,: m„ = 3M+^f2-C-ln(-/-)') , (3.7) 

where C ~ 0.577 is Euler's constant. The spatial distribution of the baryonic 
correlators is Gaussian in both p and if, and have widths ~ ATr/ey/M. Another 
simple alternative to a full solution to the Faddeev equations is to use a quark- 
diquark model as shown in Appendix C. 

Like in the mesonic case, the degeneracy of the nucleon and delta screening 
masses will be lifted by spin dependent effects. At very high temperatures we 
expect the splitting to be due to the spin-spin interaction (3.3). This gives an 
upward shift to both the nucleon and the delta. The removal of the degeneracy 

^ As discussed in I, this theory, while of no direct physical interest is easier to handle 
theoretically and perfectly amenable to lattice calculations. 

8 



is, however, only partial as the A(3/2, ±1/2) states remain degenerate with 
the iV(l/2, ±1/2) states. We get 

±1/2 ±1/2 1 / ±3/2 \ /„ ox 

m^ -ma=m^' - rria = - [m^ - rriaj (3.8) 

and the splitting is ~ e^/4M^. As already mentioned, the different polariza- 
tion's of the vector and axial vector mesons, are split so the j = ±1 states 
are pushed up, while the j = state is unaffected along with the pion and 
its scalar partner the a. At very high temperature where this effect should be 
dominant, we thus predict 

±3/2 ±1/2 ±3/2 ±1/2 /„ n\ 

m^ — m^ = m^ — rrij^ (3.9) 

of the same order of magnitude and direction as the meson splitting 

m^^ — m° = m^^ — m,r (3.10) 

Since chiral symmetry is manifest at high temperature, these relations also 
hold for the chiral partners. 

4. Gluonic Correlators 

In this section we discuss correlators of local operators containing gluon 
fields. Color singlets are easily constructed from the field strength tensor Gy^y^ 
and its covariant derivatives. We shall only consider the operators E"^ and 
5^, which are the lowest dimensional operators that couple to the two gluon 
channel, and the trace of the Polyakov loop. 



where P denotes path-ordering, which has been extensively studied by lattice 
calculations. 

At finite temperature, E'^ and B^ are distinct operators since Lorentz in- 
variance is broken. In the dimensional reduction scheme, this is manifest by 
the presence of both a (2+1) dimensional vector gluon field, A^ and a scalar 
adjoint Higgs field, if, in the lagrangian (2.5). The operator E'^ couples to the 
two-particle channel {HH) and (A^A^) with strength 1 and g'^ respectively, 
while B^ only couples to (A^A^), see fig. 1. Since g -^ at high T, we would 
thus naively expect the screening mass 2'mE in the [E'^ E'^) correlator and a 
power-law fall off in the (B"^ B^) one. Of course we know that this is oversim- 
plified. As already discussed there are lots of evidence that the dimensionally 

9 



reduced theory is confining so it makes no sense to treat the magnetic sector, 
i.e. the A^ field, in perturbation theory. We might hope, however, that the 
(HH) channel can be treated in the same way as the qq one. In that case the 
above asymptotic [i.e. T -^ oo) prediction 

{E\x\ x^)E\0, 0)) ~ Ki{mE2X^) (4.1) 

where mE2 = 2mE = -j^gT, and Ki a modified Bessel function, should be good 
at least for some moderately large values of X3. (For really large distances there 
will always be an uncontrollable contamination by the magnetic sector since 
the small coupling strength will be compensated for by a larger exponential 
factor. This is not the case in 3d QCD at high temperature where we expect 
the above result to hold with the substitution Ki —>■ Kq). It is important to 
stress that if it were not for the correlations in the transverse directions, (4.1) 
would decay as Kf{mEX^), which has the same exponential fall-off as (4.1). It 
is the pre-exponent factor that distinguishes between bound and free electric 
gluons at asymptotic temperatures. 

The binding energy is easily calculated using formulae similar to (2.7) and 
(2.8) modified for the particles having octet color charge and being scalars. A 
simple test of the self consistency of the estimate (4.1) is that the radius of 
the corresponding bound state is small enough for the magnetic effects to be 
ignored. Using the same variational estimates as in the mesonic case yields 
bound state radii ~ 1/g^^'^T, which is in between the electric and magnetic 
length scales 1/gT and 1/g'^T respectively. We consider it an open question 
whether the method based on dimensional reduction is reliable in this case. We 
should also mention that E"^ also couples to three and four particle channels, 
but these contributions are again suppressed by exponential's and/or powers 
of ^. 

We already concluded that the operator B^, is beyond the range of appli- 
cability of our method. If, however, the main effect of the infrared divergences 
in the magnetic sector is the generation of a magnetic mass rriM ~ g^T, we 
have the qualitative prediction that asymptotically the B^ screening length is 
larger than the E'^, since the magnetic mass is smaller than the electric. It is 
clearly of interest to simulate both these correlators on the lattice. 

We now turn to correlators of Polyakov loops (or Wilson lines) , 

(TrLt(f)TrL(O)) -^ iv2e-^i(^'-)/^ (4.2) 

corresponding to inserting two static quarks in a heat bath of gluons. Vi in 
( ^.2| ) is the singlet potential. In the confined phase the singlet potential is 

10 



believed to behave as a{T)x for large x = \x\, where cr(T) is the (tempera- 
ture dependent) string tension. In the deconfined high temperature phase the 
potential is usually believed to be screened, 

Vi{T,x) = 2Fq{T) + ri,o(r)e-2'"^" . (4.3) 

The exponential fall off is again governed by the two (electric) gluon channel 
just as for the E"^ correlator. As pointed out by Nadkarni []T^, also in this 
case we expect contamination from the magnetic sector at large enough dis- 
tances. However, using the same arguments as for E'^ we suggest that the 
screening length should be calculable in the dimensionally reduced theory as 
an {HH) bound state. Clearly, one would need very precise measurements to 
test whether the non-relativistic bound state interpretation of the gluonic cor- 
relators is correct. In particular it would be interesting to simulate the above 
correlator on the lattice and compare the results with the Wilson like loop at 
finite temperature. 

Recently there have been lattice simulations aimed at probing the light 
fermion distributions {ipip and ip'^ip) [|l^. Specifically, it has been observed 



that the fermion number distribution around a heavy quark given by 

(TrL(O)V'V(^)) (4.4) 

is localized around the heavy source (with a weight of minus one) at low 
temperature but spread out at high temperature. We expect that this spread 
is related to the screened two-gluon (dipole) exchange with a range of the 
order of l/2inE as discussed above. The correlation function (4.4) is consistent 
with Gauss' law in the following sense. In the presence of light fermions the 
sampled gauge configurations on the lattice have an excess or a deficiency 
of color charge (in the sense of a grand canonical colored ensemble). When 
introducing a heavy triplet source, the configurations with an anti-triplet or a 
sextet of color are sampled out so that on the average Gauss law is enforced. 
The sampled configurations together with the heavy source carry zero triality.[| 

5. Nonasymptotic Effects 

As already mentioned, our estimate of the fine structure of the screen- 
ing mass spectrum is only reliable at asymptotically high temperature. Close 
to the transition region, where the available lattice data are collected, many 
other effects are possible. Examples are Higgs interactions, effects induced 
by a dilute gas of monopoles, a residual rarefied gas of instantons and antiin- 
stantons, nonvanishing magnetic condensates... Such effects are unfortunately 

^ We thank Carleton DeTar for a discussion on this point. 
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much harder to calculate in a controlled manner, but we shall nevertheless 
attempt some estimates. 

At zero temperature, the presence of gluon condensates modifies the pertur- 
bative spin-spin interaction and gives a contribution to the spin-spin splitting 
~ (-B^) [|Tn],PD|. These effects will presumably persist even at finite temper- 



ature, but since we lack information about the condensates they are hard to 
estimate. However, the symmetry structure of these terms would be the same 
as for the perturbative spin-spin interaction (2.8). 

The effects due to Higgs exchange is rather straightforward to estimate. 
First we consider the m = case where the Higgs couple only to the scalar 
density (pcf). To lowest order, we just have to consider the exchange of a scalar 
octet particle with mass rriE- The coupling is only to the color charge, so this 
will not cause any level splitting. Rather than trying to solve the Schrodinger 
equation in a potential including the Yukawa potential from the Higgs ex- 
change, we consider the latter as a perturbation. If we furthermore neglect 
the propagation of the Higgs particle, i.e. make the local approximation l/m'^ 
for the propagator, we get the following local potential, 

Vm,,s = ^S\r) = b\^ , (5.1) 

where the last equality follows for two light quarks. This interaction will shift 
the meson screening masses by 

Am = ||^(0)p (5.2) 

and comparing with the n — p mass-shift due to the spin-spin splitting 



we see that the Higgs effect is an order of magnitude larger at a^ = 0.1 . 
Similar estimates of the level-shifts due to Higgs exchange can be made for 
the baryon and gluon correlators. 

In order to split the rho from the pi without splitting the polarization 
components of the rho, we need a coupling to 0730. A possible, and even 
plausible, scenario is that the effective lagrangian for the quarks includes an 
effective four fermi interaction of the Nambu Jona-Lassinio type. Below the 
transition temperature such an interaction is known to induce a spontaneous 
breaking of chiral symmetry (at the mean field level) and generate a constituent 
quark mass. It is possible that the same interaction could be responsible for 

12 



the large observed splitting of the pi and rho screening masses. Neglecting 
flavour (which is anyhow not correctly represented in the lattice calculations) 
we can take 



^NJL — A 



(0750) - (0730)' (5.4) 



which is chirally invariant. The corresponding spin-spin interaction between 
a quark and an antiquark is of the form ~ \S\S2 which is appropriate for 
pushing down the pi and leaving the rho untouched. The interaction is how- 
ever too singular to be treated in perturbation theory. A correct treatment, 
which would amount to resolving the Schrodinger equation with a pointlike 
interaction, will not be attempted here. 

It is also instructive to consider the limit M ^ T (heavy quark limit). 
In this case the coupling to the Higgs is entirely due to the density i0730, 
which reduces to l/2Mx^{^^-'qi(Jj)x ^^ the nonrelativistic limit, where x is a 
two-spinor and g'the momentum transfer. The resulting spin-spin interaction 
due to the Higgs is of the form 



The local part in (6.5) combines with the local part of (2.8) to give 



. (5.5) 



^^« = E l4 (1 + 4^^ ■ S^) S\r- - f,) (5.6) 

i<j 

So, as expected, for heavy quarks, we recover the conventional polarization 
independent spin-spin interaction. This result is also supported by the lattice 
results. For mesonic systems we have Hgs = (25*^ — l)e^/4M^, while for 
baryonic systems we have Hss = (25*^—3/2)6^/871/^. This spin-spin interaction 
pushes the rho up, the pion down, leaves the nucleon unchanged and pushes 
the delta up. Specifically, rrip — rrij^ = e^/M^ and m^ — -m^v = 3e^/4M^. In 
this regime we expect 

3 

rriA - rriN = -{rrip - m^) . (5.7) 

It is not ruled out, that the lattice results for the sources quoted are carried out 

with still large quark masses, in which case the present discussion is pertinent. 

Finally we give a rough estimate of the effect of instantons assuming that as 

chiral symmetry is restored, the instanton-antiinstanton medium gets rarefied 
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(single particles, molecules, small clusters). The effects would be to push down 
the nucleon but leave the delta unaffected. Since the instanton contribution is 
strongest in the quark-antiquarlc or quarlc-quarlc spin-isospin zero channel, we 
would expect a delta-nucleon splitting in the screening lengths in comparison 
to the p — IT splitting of the order of 

3 

TTiA - rriN r^ -{mp - m^) . (5.8) 

Above the transition temperature this splitting is entirely due to the nonva- 
nishing of the strange quark mass. The splitting vanishes exponentially at 
asymptotic temperatures. This splitting is comparable in magnitude and di- 
rection to the spin-spin splitting discussed above for heavy quarks. 

6. Discussion 

We have presented simple estimates of the baryonic screening lengths and 
spatial extension using arguments on dimensional reduction. These arguments 
are straightforward extensions of our previous arguments for mesons. We 
have also provided some estimates for the gluonic correlation functions, and 
suggested new correlation functions to be evaluated on the lattice. 

While the dimensional reduction scheme is only valid at asymptotic tem- 
peratures, we have also used it to estimate nonasymptotic effects caused by 
both perturbative gluons (spin-spin interaction) and nonperturbative gluons 
(instantons). It would be interesting to see how these predictions compare 
with lattice simulations for the hadronic correlation functions. We stress that 
our results are only suggestive since the calculations may be sensitive to the 
magnetic length scale. If, however, future lattice calculations will confirm, or 
at least support, our calculations, our calculations of nonasymptotic effects 
could be much improved. Since the problem is essentially that of bound states 
in 2+1 dimensional QCD, we can imagine to use many of the phenomeno logi- 
cal methods, like sum rules, bags etc., that have been developed to calculate 
bound states in 3+1 dimensions. It remains to be seen whether such a "high 
T bound state phenomenology" will become feasible. 

There are several places for improvement in our calculations. Most note- 
worthy: if our estimates of sizes and binding energy are correct, they are 
midway between the electric and magnetic scale. Thus, it is important to do a 
more complete calculation including hard thermal loops to better assess these 
effects. Since all the correlators discussed in this work are gauge invariant, the 



^We recall that the splitting is proportional to the instanton density which is proportional 
to e~'^'^^ P ^ where p is the instanton size. For temperatures of the order of Tc the density 
drops by an order of magnitude pi| , p3] . 
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tricky issue of gauge dependence related to the subdivision into hard (g'^T) 
and soft (gT) does not arise. This exercise is worth pursuing in both four 
and three dimensional QCD. Again, lattice measurements are crucial to settle 
these questions. 
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Appendix A. Two-component formulation 

In this appendix we shall show how to formulate the dimensionally reduced 
theory in a two-component formalism. We again start from the Lagrangian 
(2.1) (for simplicity we only consider the lowest positive node), but we now 
rotate to Minkowski space by the transformations 

^t^^^ = ^^t^O (^j^^^^ 

where 7° = —73. We also introduce the transverse Minkowski gamma matrices 

by 7i = iY to get 

C = V^(i7°(9o + i^di - ifM)ilj (A.2) 

After introducing the two two-spinors ipi and 1^2 via 



V- = ' (A.3) 




the Lagrangian takes the form 

C = i^.iYA^i + ^2illdi^2 - M(4^,^2 + V^2^i) (A.4) 

where the 3-d gamma matrices are defined by (73,73,73) = {cr^ , icr^ , icr"^) , and 
V'j = Vili- We leave it as an exercise to put in the gluons and a mass for 
the quark. (A.4) has the form of a pair of standard 2-1-1 dimensional Dirac 
Lagrangian's, but with an off-diagonal mass term. This means that in order to 
obtain a propagator one must diagonalize a four-dimensional matrix anyway, 
and it is more convenient to use the representation given in section 2. 

Appendix B. Quantum numbers of baryonic currents 

There are two inequivalent local sources for spin ^ baryons of opposite 
parity 
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B+ = U^C^)%^ (B.i: 



B^ = {^P'^C%^P)%^P (B.2) 

For convenience, we will throughout omit reference to color, flavour and space 
indices. Under parity {ip -^ 74^") : B^ -^ ±74!?^. B~^ and B~ refer to 
the octet | and | baryons respectively. The nucleon is an arbitrary linear 
combination of -8^25 ^•^- ^^ = -^i^ + i^B2. 

In the chiral basis the combinations B^ ± B2 are more appropriate to 
use. In this basis the nucleon can be viewed either as (RRR) — {LLL) or 
(RRL) — [LLR) or any linear combination of the two. The odd parity partners 
of the nucleon are (RRR) + [LLL) or (RRL) + {LLR) or again any linear 
combination of the two. This representation is suited to the two component 
formulation discussed in Appendix A. 

To exhibit the structure of the currents in the four-component formula- 
tion discussed in the text, we first derive the reduced forms for the di quark 
constituents. Specifically 



P^^C%^ = -0!:7i73f/!0^e-*^-^"^-0^7i73f/+0+e*'"'^"^ (B.3) 

- 0^71730+ - 0+71730- 

P^^C^ = 0^72740- e~^2.T£4 + ^T^^^^^^e.2.rx4 (3.4) 

+ 0^7274f^+0+ + 0+7274f^^0- 

The unaveraged sources for the baryons follow from these diquark structures 
and the reduced fermionic fields (2.3). The resulting expressions involve terms 
with the r-dependence e^^^'^'^'^ and e^*'^"^'^. Static sources can be constructed 
by averaging over vrT or SvrT. We choose the latter and define 

5f 2 = T r dx4 cos{37TTxi)B^2 (B-5) 

JO 

the lattice calculations have been carried out by averaging over the lowest 
Matsubara frequency vrT. Using ( [B.5| ), the reduced diquark sources (B.3) and 
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(B.4), and the fermionic fields (2.3), we obtain 

2f3lBt = -2 (0^710-) V^f/-0-+^ (0+710+) V^f/+0+ 

2/5^5+ = (0^72740-) 75l^t/-0- + (0+72740+) 75 V^t/+0+ (B.6) 

2/3iS- = (0^72740-) Vf/-0- + (0^72740+) V^f/+0+ 

2plBt = -i{4>-lik)%VUj.+t{4>llik)%VUj+ (B.7) 

Our static sources do not mix the 0+ and 0„ fields while the lattice sources 
averaged over cos(7rTa:4) do. It would be interesting to find out how (if at all) 
this affects the structure of the resulting correlation functions. 

In the case of the isobar current, there is a unique local current (■0"^C7^-0)-0. 
Its form in the four component formalism can be obtained without difficulty 
using techniques similar to the ones above, and will not be given here. We 
note, however, that the isobar carries a polarization. The isobar diquark field is 
just ■0"^C*7^-0. This means that the sources for the transverse and longitudinal 
fields are distinct at high temperature and are subject to the same remarks 
developed in the text for the vector meson fields. In particular, we expect 
the subleading spin-spin interactions to affect the longitudinal and transverse 
parts of the isobar differently. 

Appendix C. Quark-diquark model for baryons 

If we were to use the quark-diquark picture for baryons than the screening 
lengths and wavefunctions are readily obtained from the mesonic screening 
lengths and wavefunctions. In the case of a free diquark of mass M^ = 2M 
with antitriplet charge, a variational estimation shows that the screening mass 

is, 

QCD,,,: ,„„ = 3M + -^^ (C,l) 

QCD3« : m„ = 3Af + |1 (2 - C - ln(^)) (C,2) 

In the antitriple channel, the diquark mass is, 

QCZ^2+i : M. = 2M+l (^^) ' (C.3) 

gCD3+i: M^=2M+^(2-C-\n{^)\ (C.4) 
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The baryonic screening lengths in this case are, 



«.... ™.-^^.M..|(|l(^^)' (C, 



e^ / „ . , e2 



QCT„: „„ = M + M. + -^2-C-l„(^^^^^^^)J (C.6) 

In the diquark picture, the baryonic wavefunctions are similar to the mesonic 
wavefunctions. 
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Figure Captions 



Fig. 1 : Typical couplings to E"^ (a) and B"^ (b) to tlie electric and magnetic 
gluons. 



Fig. 2 : (a) Direct contribution to the sigma-sigma correlator; (b) exchange 
contribution to the sigma-sigma correlator. 
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